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The unbinding of vortex-antivortex pairs for the classical two-dimensional XY model in a magnetic 
field is studied. A single such pair is connected by a string of overturned spins, leading to linear 
confinement. We show that this system supports two phase transitions, one in which closed strings 
proliferate, and a second in which vortices unbind. The transitions are shown to be dual to one 
another, and are remarkably continuous. Possible consequences for a variety of systems are discussed. 
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Introduction. One of the most studied systems 
in condensed matter physics is the two-dimensional XY 
model, 

H/k B T= -K cos(9i - Oj) - h^cosOi. (1) 

<ij> i 

In Eq. [j], 6i represents the angle with respect to some 
reference direction of a planar spin located at lattice site 
i, J2<ij> denotes a sum over nearest neighbors, and h 
represents the effect of a magnetic field tending to align 
spins along 9i — 0. This model receives so much attention 
because it, and the closely related sine-Gordon model 
, describe a diverse range of two-dimensional classi- 
cal, and one-dimensional quantum, systems. The model 
supports topological excitations in the form of vortices, 
which for h = and high enough temperature unbind via 
a Kosterlitz-Thouless (KT) transition. 

Surprisingly little work has been done to study the 
fate of this transition for h > 0. Early renormalization 
group (RG) studies || performed on the above model, 
with the last term replaced by hJ2i C0S {P@i) with p an 
integer, indicated a range of parameters for which the 
term is irrelevant, presumably allowing a KT transition 
to be realized. For p < 4 the symmetry-breaking term 
is always relevant and the RG results are inconclusive. 
Nevertheless, for p = 1, the spins trivially exhibit long- 
range order, and it is commonly supposed that there can 
be no transition in this circumstance [Q. 

In this article, we demonstrate that vortex unbinding 
does indeed occur in the XY model with a magnetic 
field, although the nature of the transition is consider- 
ably changed from the KT behavior. In the absence of 
fluctuations, for large separations a state with a single 
vortex-antivortex pair will contain a string of reversed 
spins connecting them, leading to a linear pair poten- 
tial ||. In the absence of fluctuations, the width of this 
string is given by £o = a Q\/K/h, where ao is the under- 
lying (square) lattice constant, and the energy per unit 
length is 8V Kh (fc^T/ao) 0. Even for small K, pairs 
will be tightly bound and relatively isolated if h and E c , 
the vortex core energy, are large. As h and E c decrease, 
vortices unbind in a two-step process: a first phase transi- 
tion occurs in which the vortices remain bound in pairs, 
but the strings connecting them become unbounded in 



length. The resulting state may be thought of as one in 
which closed strings have proliferated. This is followed 
by a transition in which the strings break open, with 
the endpoints being the vortices. This sequence has a 
dual description in terms of a solid-on-solid model with 
screw dislocations 0. Interestingly, the vortex unbind- 
ing maps to a rough-flat transition (equivalently, domain 
wall proliferation) in the dual model, and string prolifer- 
ation maps to dislocation unbinding. The transitions are 
remarkable in that the free energy is perfectly analytic as 
the system passes through them, in spite of a diverging 
length scale and a clear change in the low energy dynam- 
ics across the transition. Fig. 1 illustrates our proposed 
phase diagram for this system in the E c vs. 1/h plane, 
for a fixed value of K « l/2n. 

Such deconfinement transitions have consequences for 
a variety of systems, some of which we will briefly de- 
scribe towards the end of this article. 

String Model. We begin by explaining why in a mag- 
netic field, vortices come with strings attached. For this 
purpose, we write the continuum version of Eq. |l|, 

H = J d 2 r{^K\V9(r)\ 2 - hcos9{r)} 7 (2) 

where we have chosen units in which fc^T = 1 and ao = 
1. A vortex state is a local energy minimum of Eq. ||, 
subject to the constraint that J p dr ■ V6* = 27m, n an 
integer, along any closed path P surrounding the vortex 
core. Consider a state with a single vortex in a finite size 
system, and a large path surrounding it. Because of the 
field, nearly all the spins along the path must point in 
the preferred direction; the 27rn rotation will arise as a 
localized kink. If we choose the path parallel to W9 at 
the location of the kink and choose coordinates with the 
y direction parallel to the path, |V#(f)| 2 « (d9/dy) 2 , and 
we recognize the kink as the soliton of the sine-Gordon 
model ||. Since this kink must appear for any large 
path enclosing the vortex, it is clear there is a string of 
overturned spins running from the vortex to the system 
boundary. The existence of such string defects in the 
context of the XY model in a magnetic field was noted 
in Ref. §. 

The unbinding transition for h > is driven by fluctua- 
tions of the strings. To see this, we replace the partition 



1 



function for Eq. [l] by the Villain model, in which the 
substitution 



„Ccos0 



E 



-C(e-27i-s) 2 /2 



(3) 



is made when cosines appear in the exponent. The Gaus- 
sian dependence of 9 in the Villain model allows some 
progress in evaluating the partition function without sac- 
rificing the critical properties of the XY model, since the 
two models share the same discrete periodic symmetry. 
Using the techniques of Ref. M, the 9 field may be in- 
tegrated out, and with some manipulation || the par- 
tition function for the Villain model may be written as 

Z V M = E{n(r),A(r)} eX P[-^M], with 
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In this expression n(r) and A(r) are integer degrees of 
freedom fl0|| , and derivatives should be understood as 
really meaning discrete differences; e.g., dA/dy = A(r + 
y)-A(r). 

It is useful to think of the n variables as residing at 
the centers of the square plaquettes formed by the lattice 
sites, and the A variables as residing on the nearest neigh- 
bor bonds in the y direction. Hvm may then be thought 
of as a solid-on-solid model with screw dislocations 0. 
In this interpretation, n is an integer height variable, and 
non- vanishing derivatives dn/dx, dn/dy locate domain 
walls between different heights. The important low en- 
ergy excitations involving A / contain line segments 
along the domain walls in which dn/dx + A = 0, effec- 
tively removing part but not all of a domain wall (see Fig. 
2). An endpoint of an "A" string - where dA/dy ^ 
may be identified as the core of a screw dislocation, and 
we see the second term in Hy m actually specifies its core 
energy to be l/2h. Our theory is thus effectively one in 
which the partition sum is over graphs containing open 
strings, as might be expected from the considerations dis- 
cussed above. 

Some comments are in order before we analyze Hvm- 
The model is closely related to one introduced in Ref. 

which expresses the partition sum in terms of two 
interacting vortex fields. That model may be related to 
Zv m by eliminating n using the Poisson summation for- 
mula, and representing the A sum in terms of dA/dy 
The resulting Hamiltonian is identical to the vortex 
Hamiltonian of Refs. [|g] (see Eq. 2.75b of Ref. g) with 
E c = 0. A non- vanishing core energy may be introduced 
in that representation, after which the resulting theory 
enjoys a duality: upon interchanging l/2h and E c , and 
changing K — > \/A / k 2 K 1 the partition sum is identical to 
its original form. This observation will be important in 
obtaining the phase diagram. 

RG Analysis. To make further progress we need a 
model with continuous fields but the same symmetries as 



Hvm- A standard way to proceed 11 1 is to exchange the 
integer fields n and A for continuous variables <j) and a, 
while adding terms that favor integer values of the fields. 
Taking the continuum limit, our effective Hamiltonian is 



1 i, x , \-i2 i 2 ( da 
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Heff= Jd 2 r 

- y cos{2it4>{v)) + y- y a cos{2ira{v)) + y a 



(5) 



For large E c , y = e~ Ec is the usual vortex fugacity Jl2| . 
An important observation is that the low energy config- 
urations of H e ff are highly analogous to those of Hvm, 
even for small y, y a . Domain walls in n are analogous 
to soliton configurations of <f>, and the a-field introduces 
local minima of H e f f that are analogous to open domain 
wall configurations. An endpoint of an open wall is anal- 
ogous to a dislocation core in Hvm, by matching the 
endpoint energy to the dislocation core energy l/2h, one 
may assign a value to y a such that there is a reason- 
able mapping of the low energy configurations of Hvm 
to those of H e ff. A rough calculation yields y a « l/16/i 

Our analysis proceeds from the small y, y a limit so 
that the cosine terms may be treated perturbatively. We 
wish to integrate out short length scale fluctuations in the 
partition function, which is accomplished by integration 
of </>(k), a(k) in the shell A/6 < \k x , y \ < A, where b = e l is 
the usual RG scale factor |ll[] and A = ir is the momen- 
tum cutoff. The scaling transformation takes the form 
k = k'/b, r = br', (f>'(r') = 0(r), and a(r) = a'(r')/b. 
The renormalization of the y a J d 2 r cos [27ra(r)] term un- 
der this RG transformation is unusual and leads to the 
remarkable critical properties of the deconfincmcnt tran- 
sition. Since the argument of the cosine shrinks as the 
short wavelength fluctuations are integrated out, one is 
led to expand the cosine in a Taylor series, 

oc 

y a cos(2na) - y a = £ -g- (-1)"(2™) 2 ", (6) 

n— 1 ^ ' 

with the initial values y2n(l = 0) = y a . Upon rescal- 
ing, power counting suggests y% n ~ b 2 ~ 2n . Thus we ex- 
pect all the terms to be irrelevant except the Gaussian 
one, which neither grows nor shrinks. The resulting fixed 
point Hamiltonian has the form 



H t . 
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The parameter p plays an important role: it intro- 
duces a finite energy per unit length for creating the "a" 
strings that allow dislocation pairs to populate the sys- 
tem. Thus p is effectively a (renormalized) string tension, 
and for p > dislocations are confined. The behavior of 
jff* reflects a qualitative difference when there is a finite 
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tension: for p = there is a gapless mode for every value 
of k x when k y — > 0, whereas for p > there is only a 
single gapless mode at k = |l3|]. As we will see, the 
RG flows determine which microscopic parameters allow 
a non- vanishing p. 

In the vicinity of the scaling equations to lowest 
order in y, y n become 
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The last of these equations has been evaluated at = 
0, and r\ is a number of order unity. Note that j/2 = 47r 2 p, 
and the initial values ym(l = 0) arc all given by y a . The 
function £ arises from the momentum shell integral which 
may be evaluated exactly. For small values of p, it has 
the form 

KA 2 



■KyjKp(l + £ 2 A 2 ) ' 



A remarkable feature of Eqs. [T(j is that their solutions 
are rel ated by y 2n +2{l) = V2 n e~ 21 = 2/2n£(0 2 /£o where 
£o = \J K/h is the initial value of £(/). Eqs. |l^ can then 
be rewritten in terms of a single scaling equation, 



(12) 



Fig. 3 illustrates the trajectories resulting from Eq. [LJ 
for different initial values of £o- (The initial value of 
y a may be written in terms of £o using the core energy 
matching discussed above.) For small values of £o, p{f) 
lies below the separatrix (shown as a heavy line), touches 



the p = axis at a finite value of £ 



-21' 



, and 



remains on this axis as it flows to the origin. We as- 
sociate l scr — aoe 1 with a screening length: for sepa- 
rations below this the dislocations appear to be bound, 
while for larger separations the interaction is screened 
by other dislocations, allowing for an unbinding transi- 
tion. This length scale diverges as the phase transition 
is approached. Above the separatrix, the flows end at a 
non-vanishing value of p, which we identify as a renor- 
malized string tension: the dislocations are confined if £o 
is large enough. Thus, the separatrix represents a decon- 
finement line, with dislocations unbound if £o is smaller 
than some critical value £ cr . 

Our RG analysis has demonstrated that for very large 
values of E Cl there is a phase transition in which dis- 
locations unbind. The effect of vortices on this transi- 
tion enters at second order in y, and one can show M 



that this decreases the value of £q = K/h at which the 
transition occurs as E c decreases from infinity. On the 
other hand, at E c /K = 0, K/h = 0, it may be shown 
that HyM has the character of a set of decoupled one- 
dimensional solid-on-solid models ||, which should not 
have any phase transitions. Thus, we presume the dislo- 
cation unbinding line strikes the K/h = axis at a finite 
value value of E c /K, as shown in Fig. 1. 

The dislocation confinement transition may be thought 
of as one in which domain walls go from open strings 
to proliferated closed loops. The confined phase is thus 
analogous to the rough phase of a solid-on-solid model 
{|.[l4|. For larger values of y (smaller E c ), a second tran- 
sition to the smooth phase may occur |l5|] . In the dual of 
HyM, this transition is equivalent to vortex unbinding. 
Because of the duality in the partition sum, this transi- 
tion has precisely the same character as the dislocation 
unbinding transition, and we expect for large values of 
K/h that it should occur when E c /2tt 2 K — £ 2 r ; i.e., at 
the same value one finds for dislocation unbinding. This 
is the origin of the lower transition line in Fig. 1. It 
is interesting to note that the confined phase is actually 
one in which closed strings (loops) are proliferating. In 
the dislocation language, they are the domain walls as- 
sociated with the height n; in the dual representation (in 
which the vortex degrees of freedom are explicit) they 
correspond to strings of overturned spins. 

A remarkable property of the deconfincment transition 
is that it is perfectly continuous. The "accumulation 
point" at p = 0, £ 2 = in Fig. 3 appears at the end 
of a line of fixed points, but there is no relevant direction 
(in the RG sense) leading away from it. Such relevant 
directions are the source of nonanalyticity in the free 
energy in standard second order phase transitions [fj"l|| . 
Experimentally, this implies that thermodynamic mea- 
surements cannot detect the deconfinement transition - 
one must instead directly probe the dislocations, vortices, 
or the correlation functions they affect. 

Finally, we note that measurements of the vortex dif- 
fusion constant in XY model simulations have confirmed 
that vortex unbinding indeed occurs in this system |ll| . 

Applications. The XY system in a magnetic field and 
several closely related models may be used to describe 
a large number of systems, and the transition discussed 
in this work has consequences for most of them. These 
include bilayer thin-film superconductors, for which the 
staggered conductivity would be controlled by a decon- 
finement transition; bilayer crystals, for which disloca- 
tions are analogous to vortices of the XY model, and de- 
confinement may lead to an unlocking transition of the 
layers; and double layer quantum Hall systems at to- 
tal filling factor v = 1 ||, where vortices carry physical 
charge ±e/2 and deconfinement may lead to charge frac- 
tionalization. One dimensional quantum systems (Lut- 
tinger liquids) are also often described by models closely 
related to Eqs. [j] and ^ and the transitions described 



3 



in this work have consequences under a variety of circum- 
stances. Finally, generalizations of the physics discussed 
here to layered systems, including vortices in high tem- 
perature superconductors, highly anisotropic crystals, 
and quantum fluctuations in 2+1 dimensional striped 
systems, are possible and are likely to lead to a variety 
of new phenomena. 
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Hvm accesses this transition, because the the constraints 
of duality suggest a minimum value of E c is required for 
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FIG. 1. Schematic phase diagram for XY model in a mag- 
netic field for K ~ 1/2-7T. Upper left corner represents an 
ordered phase (unbound screw dislocations in dual represen- 
tation). Middle phase contains proliferated loops in both de- 
scriptions, and lower right corner contains unbound vortices 
(flat phase in dual representation.) 



Domain Will 



































































Dislocati 
































A=l! 


all -~" 




































Dislocati 















on core 



FIG. 2. A low energy configuration involving A 7^ 0. A re- 
gion of n — 1 (hatched squares) is embedded in a surrounding 
n = region (white squares); heavy line represents a domain 
wall. The line segment with ^4 = 1 cancels the domain wall 
energy for part of its length, leaving an open domain wall. 
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FIG. 3. RG flows for Eq. Left vertical axis is a fixed 
line, and heavy line is a separatrix between flows that reach 
p = for finite I and flow to the origin as I — > oo, and those 
that have p > at the end of their flow. Finite values of p in 
H, indicate confinement. 



